Broad Ha wings in Nebulae around Evolved Stars and in Young 

Planetary Nebulae 

A. Arrieta 1 AND S. Torres-Peimbert 

Institute) de Astronomia (UNAM), Ap. Postal 70-264, Mexico D.F. 04510 Mexico 

q 1 anabeKaastroscu.unam.mx 

O ' 
(N 

X> ■ ABSTRACT 

^ ! Eleven objects that have been reported as proto-planetary nebula or as young 

planetary nebulae that show very extended Ha wings are presented. The exten- 
sion of these wings is larger than 800 km s _1 . Data for two symbiotic stars that 
show this same characteristic is also presented. Raman scattering is the mecha- 
^ ! nism that best explains the wings in 10 of the PNe and in the 2 symbiotic stars. 

In the PN IRAS 20462+3416 the wing profile can be explained by very intense 
\ stellar wind. 

O 

Subject headings: ISM:planetary nebulae:-binaries:symbiotic-scattering- stars:mass- 



O 



loss-stars: AGB and post-AGB - Line: profiles 



1. Introduction 



As part of a study of proto-planetary nebulae and young planetary nebulae we have 
observed 59 objects that have been selected for fulfilling one or more of the following criteria: 
post-AGB stars with evidence of high mass loss rate, or PNe with low degree of ionization 
or with H2, CO or OH molecular emission; all the objects observed have IRAS colors typical 
of planetary nebulae. In the sample, 13 objects were outstanding for exhibiting extremely 
broad Ha lines, of these, 10 have not been reported previously as having this characteristic. 
The Ha profiles selected for being extremely broadened are presented in this work. 

Very wide Ha emission lines have been reported previously in other objects in the early 
stages of planetary nebulae phase; these include seven post-AGB stars, five young PNe and 
some symbiotic stars (Van de Steene, Wood & van Hoof 2000; Lee 2000; Lee & Hyung 
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2000; Miranda, Torrelles & Eiroa 1996; Balick 1989; Lopez & Meaburn 1983; Van Winckel, 
Duerbeck & Schwarz 1993; Wallerstein 1978). Lee & Hyung (2000) advanced the proposal 
that the broad wings are produced by Rayleigh-Raman scattering which involves atomic 
hydrogen, whereby Lya photons with a velocity width of a few 10 2 kms^ 1 are converted 
to optical photons and fill the Ha broad region. In the present work this mechanism is 
investigated further and alternate possibilities for Ha broadening are examined. 

2. Observations and reductions 

The observations were carried out at the Observatorio Astronomico Nacional in San 
Pedro Martir, Baja California with the 2.1-m telescope and the REOSC echelle spectrograph 
(R ~ 18,000 at 5,000 A) and a 1024x1024 Tektronix detector that yields a spectral resolution 
of 10.6 km s _1 and a spatial resolution of 0.99" per pixel. The 3600 to 6800A range was 
covered in 29 orders. For extended objects the slit was centered on the nucleus. The 
observing log of those objects with wide Ha wings is presented in Table I. 

The data reduction was carried out with IRAF routines using 'echelle' and 'ccdred' 
tasks. Calibrated spectra for most of the objects were obtained. The orders were extracted 
using a 10 pixel window. A HeAr lamp was used for wavelength calibration and the standard 
stars by Hamuy et al. (1992) were observed for the flux calibration. 

In several cases the objects exhibit extremely wide Ha lines, with full widths at zero 
intensity, FWZI, larger than 800 kms -1 . In order to ascertain whether the wide wings are 
real or due to instrumental effects the "number of counts at maximum" vs. FWZI both for 
the brightest unsaturated lines of the comparison lamp spectra and for Ha emission lines 
in our sample were compared. Objects with with FWZI > 1000 kins" 1 and those with 
significantly broader wings in Ha than what can be considered instrumental were selected 
for this study; namely, those with FWZI 3<r above the observed ones for comparable line 
intensities in the HeAr lamp. 

3. Ha wings 

The Ha profiles for the 13 selected objects are shown on Figure 1. To better visualize the 
line profiles, an enlargement along the intensity to display the wings over the full Ha profile 
is presented. In all cases the FWZI is significantly larger than expected from a gaussian fit 
to the core of the profile. No correlation between the line intensity and the width of the 
wings was found. 
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The main spectral characteristics of each object are presented in Table ??. It includes 
the ions present in emission and absorption in our spectra and the stellar spectral type (as 
taken from SIMBAD database). For the 13 objects Ha FWZI and H/3 FWZI where it is 
significantly larger than expected from a gaussian fit to the core of the profile are presented 
in Table 3. In this table the most common classification in the literature for each object, 
nebular morphology, possible binarity and type of Ha line profile are included. 



4. Possible mechanisms for wing broadening 

4.1. Rotating disks 

For this case maximum velocity would correspond to the circular velocities at the surface 
of the exciting star, 

v m ax = 437 v / (m/M )/(r/ J R ) kms" 1 , 

which for young PNe central stars of m = 0.6 — 0.83 M Q , and r = 1 — 18 R Q would correspond 
to v max < 400 kms -1 , although for symbiotic stars, where the hot component is small r < 0.1 
R , larger velocities could be achieved. Nevertheless, for both PNe and symbiotic stars, v max 
is not large enough to explain the observed broadening which reaches up to 5000 kms -1 in 
some cases. Thus the possibility that in PNe emission from a rotating disk would explain 
by itself the extended wings can be ruled out. 



4.2. Stellar winds 

The possibility for the line wings to be formed in the region dominated by stellar winds 
was tested by examining the existing IUE spectra for each object looking for evidence of 
P Cyg profiles in resonant lines of ions like C IV, C III, He II, Al III, Mg II and Si IV. 

Only one object (IRAS 20462+3416) shows P Cyg profile evidence both in our optical 
spectra and in the UV. Furthermore, IRAS 20462+3416 shows an anomalous broadened 
emission in Ha, that is +110 kms -1 from the narrow emission line. This profile shows 
significant deviations from the AA~ 2 profile that can be fit to the other broad wings (see 
Figure 1). The presence of P Cyg profiles both in optical and UV spectra, and the Ha wing 
profile suggests that stellar wind is the mechanism that is producing these features. In a 
separate paper (Arrieta, Torres- Peimbert, Georgiev & Koenigsberger 2003) it is shown that 
indeed the Ha profile of this object can be explained by a very strong stellar wind. The Ha 
profiles for the rest of the objects cannot be explained with stellar winds. 
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4.3. Electron scattering 

This mechanism has been proposed, and has been intensively studied, as the line broad- 
ening mechanism in QSOs (e.g. Mathis 1970; Shields & McKee 1981; Lee 1999) and in WRs 
(e.g. Hillier 1991). Given that the cross section of electron scattering is independent of 
wavelength it is to be expected that other intense emission lines formed in the same region 
as Ha are similarly broad. 

Extended wings in other bright emission lines of our optical spectra were searched 
for. No broad extensions in forbidden lines in any of the objects were found. The objects 
IRC+10420, M 1-92, HM Sge, M 3-60, and Z And showed extended wings in H/3. In all cases 
the width of B./3 is significantly smaller than that of Ha (see Table 3). It is possible that in 
the other objects the low signal-to-noise ratio masks the low intensity extended wings. For 
objects with broad H/3 wings no correlations were found with the presence of the H[3 wings 
and the width of Ha, nor with the Ha/H/5 line ratios. 

The necessary conditions for the Ha wings to be produced by electron scattering were 
investigated. Given that no significant broadening was found in the forbidden lines, it is 
required for the forbidden line region to be exterior to the electron scattering region. A two- 
region geometry was assumed: (a) a recombination line emitting region where a substantial 
fraction of Ha is produced, dense enough to suppress the forbidden lines, and where most of 
the scattering takes place (N e > 10 6 , T e ~ 10 4 ) (b) and an external low density forbidden 
line emitting region (N e = 10 4 — 10 6 cm -3 , T e ~ 10 4 ). It is also possible to consider that 
region (a) consists of a point source where most of the Ha is produced surrounded by a dense 
scattering region. 

To derive the temperature in the electron scattering region from the observed Ha line 
profiles the treatment proposed by Mathis (1970) was followed. The main parameter is the 
width of the wings, defined as Y(2) = F obs (Av M ) / F obs (2Av M ), which essentially determines 
the temperature of the electron scattering region, T ES . The optical depth of the scattering 
region, tes, can be derived from the ratio of the strength of the unscattered core of the line 
relative to its total strength w = F {unscattered) / F (total) , which determines the optical 
depth of the scattering region, tes (Osterbrock 1980). For the simple case of a spherically 
symmetric scattering region of uniform density and temperature for the case of M2-9 Y(2) = 
2.05, and w = 0.83 can be obtained. From these values a temperature, T e ~ 10 8 K, and 
optical depth, t E s ~ 0.24 can be found, the latter coupled to a limiting size of r E s < 10 12 
cm, given by an unresolved central core of 2", yields a density N e > 10 12 cm~ 3 for the 
electron-scattering region. From the values derived for M 2-9, which is representative of 
the whole sample, it was considered unlikely for this configuration to be the general case 
surrounding an AGB star and thus to account for the extremely wide Ha lines observed in 
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young PNe. 

4.4. Raman Scattering 

Raman scattering describes the absorption of a photon, followed by the immediate 
re-emission of another photon at different wavelength where the intermediate state does 
not correspond to a true bound state of the atom. Initially this process was suggested by 
Schmid (1989) as the mechanism of production of the broad emission features at AA6830 and 
7088A found in 50% of symbiotic stars; he suggested that the emission features are due to 
Raman scattering of the O VI resonance doublet AA 1032 and 1038A by neutral hydrogen. 
Nussbaumer, Schmid & Vogel (1989) proposed a list of uv lines of different ions (including 
S III, He II, O I, O VI and C II), with wavelengths close to Ly/3, as candidates for Raman 
scattering on the ground state of neutral hydrogen; these Raman scattered lines would be 
expected to be found in the 6000 - 7000A , range. They noted the possibility that Raman 
scattering may also hold a clue to the broad line wings of Ha, occasionally observed in 
Seyfert galaxies and symbiotic stars. 

Several of the emission features at different wavelengths that can be attributed to Raman 
scattering have been identified in symbiotic stars and in planetary nebulae. In symbiotic stars 
the identifications include, among others: AA6830, 7088A from O VI AA1032A and AA7021, 
7052A from C II AA1036, 1037A in V1016 Cyg (Schmid 1989; Schild & Schmid 1996) as well 
as A4850A from He II (2-8) A972.lA, A433lA from He II (2-10)A949.3A, and A4975A from 
C III A977.0A in RR Tel (van Groningen 1993). In planetary nebulae, they include A4850A 
in NGC 7027 (Pequignot et al. 1997), and A6545A from He II (2-6) A1025A in M2-9 (Lee et 
al. 2001). 

Lee & Hyung (2000) proposed that the broad Ha wings of the PN IC 4997 are formed 
through Raman scattering that involves atomic hydrogen and, by which Lj/3 photons with 
a velocity width of a few 10 2 kms -1 are converted to optical photons and fill the Ha broad 
wing region. Their model fits the observations on the blue wing from v = 500 to 1500 
kms -1 , it requires relatively strong incident Ly(3 flux from an unresolved core of high density 
(N ~ 10 9 — 10 10 cm -3 ) and a column density for the scattering region of n#o = 2 — 4 x 10 20 
cm~ 2 . Lee (2000) further proposed that the Ha wings seen in symbiotic stars can be fit to 
Raman scattered profiles. For the optically thin case, where almost all the Ly(3 photons are 
scattered not more than once and assuming a flat incident Lj/3 profile, to the first order the 
wing profile can be approximated by a curve proportional to f{Av) = Av ~ 2 . Lee adjusted 
this curve to Ha profile observations from 200 to 1000 km s -1 from the line center for 16 
symbiotic stars. 
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Figure 1 shows the fit of curves proportional to Av " 2 to the wings of the 13 objects 
under consideration in this paper. In general the fits are satisfactory for 12 objects, excepting 
IRAS 20462+3416 where it does not at all fit the observed profile; the curves match the 
observations starting at 200 - 300 kms -1 from the center of the line, to the outermost 
regions where the signal becomes too faint to be significant, this extreme velocity reaches 
up to 500 - 1500 kms -1 . The differences in some objects can be explained in terms of non 
uniform motions in the ionized regions. The kinematics of these cases will be analyzed in 
later studies. In particular, CRL 618 and HM Sge show wider cores than the expected 
profile, while IRC+ 10420 shows a narrower core; Ml-92 shows significant differences in the 
300 - 800 km s -1 region, which have been interpreted as a jet of material along the line of 
sight (Arrieta, Torres- Peimbert, & Georgiev 2000). The case of IRAS 20462+3416 is clearly 
different, because it does not fit at all such a f(Av) = Av~ 2 profile, neither in the core, nor in 
the wings. As mentioned in the previous section the interpretation of this profile corresponds 
to a very large mass loss rate (Arrieta, Torres-Peimbert, Georgiev & Koenigsberger 2003) . 

The width required for Ly/3 emission to be responsible for the Ha wings has been 
investigated in more detail. Since the width of the scattered Ha is proportional to the initial 
width of the Ly(3 line, then, Av^yp = Ai># Q /6.4 is expected. An estimate of the Ly/3 width 
can be given by other UV emission lines widths. A search in the UV high resolution spectra 
of MAST (Multimission Archive at Space Telescope) was carried out for the objects under 
consideration where data for Z And, HM Sge and IC 4997 were found. In the cases of Z And 
and HM Sge the Si III] A1892A line has width of ~ 500 km s -1 , consistent with the Av Ly p 
width (~ 625 and ~ 468 km s -1 ) required to produce the observed broadened Ha lines (as 
given in Table 3). In the case of IC 4997 the width of the Si III] A1892A line is of only ~ 160 
km s _1 . Although Lee et al. (2001) consider that since the Ly/5 underlying continuum 
emission is not likely to show a flat profile, it may be more plausible that the broad Ha 
wings are formed by Raman scattering of the continuum photons around Ly/3. 

Evidence of additional Raman scattered features in our optical spectra was searched for. 
The A6830A feature was found to be present in Z And and weakly in HM Sge with a FWZI 
~1500 kms -1 (this feature has been studied previously by Birriel et al. 1998; Schmid et al. 
2000). Features around this wavelength were found in M 2-9, IRAS 17395-0841, IRC+10420, 
M 3-60, IC 4997 and possibly in M 1-92 with FWZI around a few 10 2 km s" 1 and not a few 
10 3 km s _1 as is the case in symbiotic stars. The A6545A feature in M 2-9 and M 1-92 as 
well as A4851A in IC 4997 with FWZI around few 10 2 of km s were also found. In that 
same sense, the H/3 wings found in IRC+10420, M 1-91, HM Sge, M3-60, Hb 12 and Z And 
(see Table 3) can be explained by Raman scattering of the Ly7 line by the neutral hydrogen 
component. 
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In order to obtain additional indicators of the presence of neutral or molecular material 
along the line of sight a search in the literature was carried out for data on 21-cm measure- 
ments and the column density of atomic hydrogen. Only in the symbiotic star HM Sge and 
the planetary nebula IC 4997 21-cm column densities have been measured; in both cases 
values of n#o ~ 4 x 10 20 cm~ 2 were determined. Also, following Dinerstein et al. (1995), 
evidence of nebular Na I doublet AA5889, 5895A from our optical spectra was looked for, 
where Na I nebular components in 11 of our 13 objects were found; those with bright opti- 
cal continuum and favorable radial velocities, (objects with radial velocities well separated 
from those expected for the interstellar material) which allow the separation of the nebular 
component and the narrow interstellar one. In addition a search in the literature for CO 
rotational transitions in radio wavelengths, roto- vibrational ones in the near infrared and H2 
roto-vibrational transitions was performed. The results are given in Table 4 where evidence 
of a neutral component at least in one of the indicators is listed. 



5. Discussion and summary 

We confirm that Raman scattering is the most probable mechanism for the formation 
of the Ha broad wings in 12 of the 13 objects under consideration. This is supported by a 
fit of the profile to a Av ~ 2 law, by the presence of Raman produced features in emission, 
and by additional indicators of the presence of a significant neutral hydrogen component. 

Evidence of neutral components in other 29 objects observed with Ha in emission was 
searched for, in order to compare with the sample presented here. However there was not 
enough information in the literature to derive significant statistical data between the objects 
that show extreme wing widening and those that do not show it. Only in the case of the 
hydrogen molecule it was found that those objects with broad Ha wings have H 2 , while those 
that are not broadened, do not. 

Most of the objects with broad lines seem to have common characteristics: bipolar 
morphology, and composite emission line profiles, furthermore, they are compact objects in 
the process of forming a planetary nebula and very probably they are sites of wind interaction. 

More detailed studies of the conditions for Raman scattering to take place and thus for 
the presence of surrounding neutral hydrogen could give us information about the processes 
that take place in giving rise to bipolar planetary nebulae. 

We are grateful to M. Bautista, L. Georgiev y E. Villaver for fruitful discussions and 
to H. M. Schmid for many important suggestions to the original manuscript. Support from 
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Fig. 1. — Ha profiles for the 13 selected objects of extreme broadening. In 12 cases a section 
of the spectrum between -2000 and +2000 kms _1 is shown, and in one case -1000 and +1000 
kms" 1 . Superposed to the observations is the 1/v 2 function that best adjusts to the wing 
profiles within ± 100 kms _1 of the core. 



Table 1: Observing log of objects with extended wings 







date 


exposure 


Object 


PNG 


observed 


time (sec) 


CRL 618 


166.4-06.5 


5 Sep 95 


1800 


M 2-9 


010.8+18.0 


23 Apr 96 


600 


IRAS 17395-0841 


017.0+11.1 


10 Jun 97 


1500 


Vy 2-2 


045.4-02.7 


3 Jun 97 


900 


IRC+ 10420 




27 Sep 95 


900 


M 1-91 


061.3+03.6 


11 Jun 97 


2700 


M 1-92 




10 Jun 97 


300 


HM Sge 




23 Apr 96 


600 


M 3-60 




27 Sep 95 


900 


IC 4997 


058.3-10.9 


9 Jun 97 


180 


IRAS 20462+3416 




6 Sep 97 


1600 


Hb 12 


111.8-02.8 


1 Sep 97 


600 


Z And 




25 Sep 95 


900 



Table 2. Nebular and stellar spectral characteristics of the objects with very extended Ha 

Wings. 



Object 



Emission lines 



Absorption lines Sp. type 



CRL 618 
M 2-9 



IRAS 17395-0841 
Vy 2-2 

IRC 10420 

Ml-91 

M 1-92 



HM Sge 



M 3-60 
IC 4997 



Hl,[FeII],[N II], [O I],[0 III],[S II] 
H I,He I,He II, [Ar III],[Cr II], [Fe I], 
[Fe II] , [Fe III] , [Fe IV] , [N II] , [Ni II] , 
[Ni III],[Ni IV], [O I],[0 II], [O III], 
[S I],[S II], Si II 

H I,He I,[N II], [Ne III],[0 I],[0 III],[S II] 
H I,He II, [Ar III],[N II], [Ne III],[0 I],[0 II], 
[O III],[S II], [S III] 

H I,Cr II,Fe I,Fe II,Mg I,Sc II,Sr II,Ti II 

H I,He I,[Fe II], [N II], [O I],[0 III], 
[S II], [S III] 

H I,He I,[Ca II], [Cr II], [Fe II], [N I], 

[N II], [O I],[0 III],[S II],Ti II 

H I,He I,He II, [Ar IV], C II,C III,[Fc II], 

[Fe III],[Fe IV], [Fe V],[Fe VII], [K IV],Mg I, 

[N II], N III,Ne II, [Ne III],0 1,0 11,0 III, 

[S II], [S III], Si II 

H I,FeII,[N II], [O I],[S II] 

H I,He I,[Fe II], [N II], [O I],[0 III],[S II] 



IRAS 20462+3416 H I,[N II], [S II], Si II 



Hb 12 



Z And 



H I,He I,He II, Ar IV,C II,C1 II,C1 III,[Fe II], 
[Fe III],Mg I,[N II] ,N III,[Ne III],[0 I],[0 II], 
[O III],[S II], [S III], Si II,Si III 
H I, He I, He II, [Cr I], Fe II, [Fe VII], 
[K V], [Mnll], [O I], Si II 



non detectable 
non detectable 



non detectable 
non detectable 

H I,Ca II,Cr II,Fe I, 
Fe II,Na I,Si II,Ti II 
non detectable 

H I,Ca II,Fc II, 

He I,N I 

non detectable 



BO 
Be 



Of 



F8 la 



Be 



B0.5IV 



M 



non detectable B 

non detectable 

H I,He I,He II,A1 III, B 

Ca II,N II,Ne 1,0 II, 

S II,Si II 

non detectable 



TiO 



M6.5 
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Table 3: Characteristics of the objects with very wide Ha wings (classification, morphology, 
binarity, profile type and widths of Ha and H/3). 



Object 


Classif. 


Morphology 


Binarity 


Profile 


Ha 
FWZI 

(kms -1 ) 


H/3 
FWZI 

(kms -1 ) 


CRL 618 


Proto-PN 


bipolar 




double 


2300 




M 2-9 


YPN 


bipolar 


probable (1) 


double 


5000 




IRAS 17395-0841 


Proto-PN 


non resolved 




simple 


800 




Vy 2-2 


YPN 


bipolar 




simple 


1400 




IRC+ 10420 


OH/IR 


bipolar 


no (2) 


comp. 


2600 


1750 


M 1-91 


YPN 


bipolar 


probable (3) 


double 


1100 




M 1-92 


YPN 


bipolar 


yes (4) 


double 


2900 


600 


HM Sge 


symb. 


bipolar 


yes (5) 


double 


3000 


1300 


M 3-60 


YPN 


non resolved 




double 


2400 


800 


IC 4997 


YPN 


bipolar 




double 


5100 




IRAS 20462+3416 


YPN 


oblate 




PCyg 


2200 




Hb 12 


YPN 


bipolar 




simple 


1800 




Z And 


symb. 


non resolved 


yes (6) 


double 


4000 


1000 



References. — (1) Schwarz et al. 1997; (2) Hrivnak et al. 1989; (3) Rodriguez et al. 2001; (4) Feibelman 
& Bruhweiler 1990; (5) Taranova & Yudin 1983; (6) Mikolajewska & Kcnyon 1996 
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Table 4: Evidence of neutral components along the line of sight of those objects with wide Ha 
lines. 

II I 



Object 


(10 20 cm- 2 ) 


Nal 


H 2 


CO 


CRL 618 




yes 


yes (3) 


yes (5) 


M 2-9 




yes 


yes (3) 


yes (5) 


IRAS 17395-0841 




yes 






Vy 2-2 




yes (2) 


yes (3) 


no (5) 


IRC+10420 




yes 




yes (6) 


M 1-91 




no 


yes (4) 


no (7) 


M 1-92 




yes 


yes (3) 


yes (8) 


HM Sge 


4(1) 


no 




yes (9) 


M 3-60 




yes 






IC 4997 


3.8 (2) 


yes (2) 


no (4) 


no (5) 


IRAS 20462+3416 




yes 






Hb 12 




yes 


yes (3) 


no (2) 


Z And 




yes 







References. — (1) Leahy et al. 1990; (2) Dinerstein et al. 1995; (3) Hora et al. 1999; (4) Kastner et al. 
1996; (5) Huggins & Healy 1989; (6) Knapp & Morris 1985; (7) Josselin et al. 2000; (8) Alcolea et al. 2000; 
(9) Mueller & Nussbaumer 1985. 



